Feedback processes in cellulose thermal 
decomposition. Implications for fire-retarding 
strategies and treatments. 

R. Ball* 

Department of Theoretical Physics, Australian National University 

Canberra ACT 0200 Australia 

A. C. McIntosh a 
J. Brindley b 

a Department of Fuel and Energy, b Department of Mathematics 
University of Leeds, Leeds LS2 9JT U.K. 

February 2, 2008 



Abstract 

A simple dynamical system that models the competitive thermokinetics and 
chemistry of cellulose decomposition is examined, with reference to evidence from 
experimental studies indicating that char formation is a low activation energy 
exothermal process and volatilization is a high activation energy endothermal pro- 
cess. The thermohydrolysis chemistry at the core of the primary competition is 
described. Essentially, the competition is between two nucleophiles, a molecule of 
water and an -OH group on Cq of an end glucosyl cation, to form either a reducing 
chain fragment with the propensity to undergo the bond-forming reactions that ul- 
timately form char or a levoglucosan-end-fragment that depolymerizes to volatile 
products. The results of this analysis suggest that promotion of char formation 
under thermal stress can actually increase the production of flammable volatiles. 
Thus we would like to convey an important safety message in this paper: in some 
situations where heat and mass transfer is restricted in cellulosic materials, such 
as furnishings, insulation, and stockpiles, the use of char-promoting treatments for 
fire retardation may have the effect of increasing the risk of flaming combustion. 

Keywords: cellulose thermal decomposition, competitive thermokinetics, dynamical model, 
char formation, fire-retarding treatments. PACS 82.40.Py 89.90 
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1 Introduction 



Thermokinetic and chemical feedback processes that determine the course and outcome 
of cellulose thermal decomposition are of particular interest in the development of appro- 
priate flame-inhibiting treatments of manufactured goods such as bedding, furnishings, 
and insulation, or stockpiled raw cellulosic material such as bagasse or hay. The fuel 
that ignites in the flaming combustion of cellulose substrates consists mostly of volatile, 
small-molecule substances that are supplied continuously by thermal degradation of the 
solid cellulose. In designing more efficient and effective fire-retardants for cellulose - 
and, more generally, in the pre-emptive control of cellulose combustion — it is thus im- 
portant to understand both the dynamics and the chemistry of decomposition. In this 
work we investigate a phenomenonological model for the thermal decomposition of cellu- 
lose substrates that is based on the following two key aspects of the known chemistry and 
thermokinetics: (1) involvement of water in promoting char formation, and (2) reaction 
enthalpy feedback to promote volatilization. 

A great many thermogravimetry (TG) and differential scanning calorimetry (DSC) ex- 
periments, complemented by qualitative and quantitative product analysis, have affirmed 
that cellulose thermal decomposition is largely a competitive process [p], ||, |3|, £|, || || [5J. 
In the big picture that has taken shape from these results and from associated modelling 
studies (for example, see || ||, [10], [□], [12|, [U| and the reviews in |TJ[ and fl5|| ) formation 
of volatiles and char are seen to be reciprocally linked to some extent. 

Some details have also emerged of the chemistry involved in these competitive pathways 
]TE| , |Hj , [THf , although a complete picture of the complex reaction network remains a 
major experimental challenge. What is clear though is that the involvement of water 
in hydrolysis reactions is crucial in determining the rates and outcomes of the com- 
petitive pathways. A thorough survey and compilation of the known chemistry and 
thermokinetics of cellulose thermal decomposition was carried out in fll9f , on the basis 
of which a dynamical model was built up that includes the role of water in the reac- 
tion chemistry. In this paper we concentrate on the kinetic and thermal aspects of the 
decomposition process, emphasizing the role played by thermal feedback in regulating 
the decomposition pathways. We demonstrate the potential for thermal instabilities 
by considering the qualitative behaviour of a very simple thermokinetic model, which 
nevertheless contains the important rate processes in cellulose thermal decomposition: 
competitive temperature-dependent reaction rates, reaction enthalpy feedback, and heat 
losses. 
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2 Basic chemistry and phenomenology of cellulose 
decomposition 

The chemistry of the primary competitive reactions is illustrated in Fig |I|. 

[Figure 1 about here.] 

The initial step is believed to be heterolytic thermal scission of glycosidic linkages at 
random chain locations in amorphous regions of the cellulose. In a dry environment the 
positively charged end (A in the figure) is rapidly cyclized to a levoglucosan end (B in 
the figure), with the -OH group on Cq of the unit as nucleophile. Thermolysis at the 
next glycosidic linkage of this fragment releases the volatile levoglucosan. When water is 
present it can compete as a nucleophile for the positively charged centre on A to produce 
a reactive reducing end (C in the figure). It is this species which is believed to undergo the 
subsequent dehydration, decarboxylation, and cross-linking reactions that produce the 
char. (The negatively charged fragment of the thermal scission rapidly picks up a positive 
hydrogen ion from the -OH nucleophile group or from water to form a non-reducing end 
D. In general, non-reducing ends of carbohydrates are relatively unreactive.) 

Empirical methods for producing good yields of char by heating damp wood in an oxygen- 
deficient atmosphere have been used for centuries. Systematic investigations of this 
effect using DSC experiments were carried out by Mok et al who observed that 
either a high concentration of vapour products or added water increased char yield and 
decreased the temperature of onset of decomposition. Subsequent analysis of these data 



20] led to the suggestion that the water produced in thermal dehydration reactions feeds 



back to hydrolyse the unreacted cellulose. In |2lJ] and |22j evidence was also presented 



from pyrolysis experiments that product water accelerated the formation of char. Some 
investigations of the role of the hydrolysis chemistry of cellulose under thermal stress in 
producing this effect were reported in [|17| and Q . 



It is also known that covered stockpiles of cellulosic materials such as bagasse, hay, and 
cotton undergo spontaneous combustion much more readily when critically damp than 
in dry conditions. Dynamical models that reproduce qualitatively the behaviour that 
is typical in damp combustion have been investigated by Gray |24| and others p5| , [26| , 
|27|, [28|, E9j. In these systems, which model the combustion subsequent to decomposition, 
critical, unstable, or oscillatory thermal behaviour can occur due to the coupling of two or 
more nonlinear temperature-dependent reaction rates with restricted linear heat losses. 
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3 Alternative pathways are thermally connected 

The studies reviewed briefly in the previous section suggest the rather counter-intuitive 
notion that water can increase the temperature of thermally decomposing cellulose. To 
guide our thinking on this, we need to put together three pieces of information: 



1. The volatilization rate is more temperature-sensitive, i.e., has higher activation 
energy, than the charring rate. A typical experimental result that leads to this 
inference was reported in |J. It was found that as the pyrolysis temperature of 
pine wood, cellulose, hemicellulose, and lignin was increased, the percentage mass 
of solid char decreased, while that of gas and oil (volatile or tar) products increased. 

2. Estimates of the magnitude of the reaction enthalpies vary, but there is agreement 
that the charring process is highly exothermal whereas volatilization is endother- 
mal, at the very least to the extent of the enthalpy of vapourization of levoglucosan 
(see for example [^1], [Jt], [H], |32|, |33|). 



3. As described in |]19fl , water of dehydration can participate in hydrolysis that effec- 
tively blocks the cyclization to levoglucosan. 



Putting these together we can draw the schematic diagram of the competitive decompo- 
sition shown in Fig [| 

[Figure 2 about here.] 

As the substrate A is heated the low activation energy dehydration reactions are initiated. 
The water produced by this process can feed back into hydrolysis reactions that inhibit 
volatilization and enhance char formation. However, the charring reactions are bond- 
forming and heat-releasing. Under circumstances where there is no efficient linear heat 
sink, the alternative, high activation energy, volatilization reactions begin to take over. 
The heat requirement of volatilization means that it can self-damp, thus opening the 
reaction field again to the charring pathway. It is clear that there is potential for the 
system to attain a stable steady state or limit cycle. In reality, of course, the situation 
is more complex — and probably more dangerous, because accumulation of hot volatiles 
in what is often a highly oxidizing environment is likely to lead to spontaneous flaming 
ignition. 
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4 A simple proxy for competitive thermal decompo- 
sition 

To strengthen our feel for the qualitative aspects of the system described above, we will 
study the thermokinetics of the following simple competitive reaction scheme: 

(1) 
(2) 
(3) 

This may be viewed as a model for the the thermal decomposition of cellulose, with 
A, W, Y, Z, and A representing the substrate, water, dehydrated product or char, 
volatiles, and heat respectively. Although a gross simplification of the real process, 
it nevertheless contains the important overall characteristics of temperature-dependent 
reaction rates, reactive heat generation/consumption, and feedback of W. Making the 
reasonable approximations that (1) the rate of supply R w of W from reaction |1] is constant 
(the activation energy for this reaction is believed to be very low) and thermally neutral, 
(2) A is drawn from an "infinite pool" , and (3) the system is spatially homogeneous with 
linear heat dissipation, we can write the following dynamical equations to describe the 
evolution of the vapour-phase components and the temperature: 

^ = -e-^w +r-fw (4) 
dr 

2± = ve -n/« _ fz ( 5 ) 

C^L = e~ l / u w - ave~^ u + £(u a -u). (6) 
dr 

The dimensionless groups and the quantities from which they are comprised are explained 
in the Appendix. Briefly, w is the concentration of reactant water, z is the concentration 
of volatile product, u is the dynamical temperature. Water is supplied by reaction [1] at 
a temperature-independent, thermally neutral rate r, and is removed with the vapour- 
phase outflow at rate / and by the hydrolysis process || that ultimately leads to char 
formation. Volatiles are produced in process [3] and are also removed in the vapour- 
phase outflow. Process is a heat source, process |3| is a heat sink, and a linear heat 
exchange with rate coefficient £ is assumed to occur between the system and its immediate 
environment at constant ambient temperature u a . 
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Eqs |]-|6] are based on a CSTR (continuous stirred tank reactor) paradigm with two 
temperature-dependent reaction rates. Comparable dynamical models that include two 
or more temperature-dependent reaction rates have been analysed in the literature from 
several points of view: fundamental studies of the dynamical and bifurcation behaviour 
may be found in [34] and |35||, applications to the wet combustion problem were cited 



above, and their use in the the stabilization of thermal runaway reactions was described in 



|36| . In all of these problems it was found that the coupling of two non-identical Arrhenius 
terms could give rise to steady-state and oscillatory behaviour of greater complexity than 
in corresponding systems with a single Arrhenius term. Concomitantly, it has been shown 
that an endothermal reaction rate (the second term on the right-hand side of Eq. ||) can 
stabilize the system temperature In this system the effects of the high activation 
energy endothermal reaction are modified by the rate r at which water is fed into the 
exothermal reaction. 

Some steady state solutions of Eqs |]-|6], computed numerically, are shown in Fig |3], where 
the dependence on u a is plotted for three values of r. Limit cycles, where they exist, 
have also been computed and plotted as the amplitude maxima against u a . (Notes: (1) 
Eq. |5] uncouples from Eqs f| and |6|. (2) With an assumed activation energy Ei of ~50 
kJ/mol a dimensionless temperature u of ~0.1 corresponds to ~600 K.) 

[Figure 3 about here.] 

In (a), where the supply of reactant is limiting, the temperature rise is monotonic as 
u a is increased. When consumption of w begins at u a ~0.08 the rate increases due to 
thermal feedback, but when production of z cuts in at a higher temperature the rate 
levels out again. 

In (b) the situation is impressively different, because the steady-state solutions change 
stability at two Hopf bifurcations. Here a faster rate of supply induces oscillatory be- 
haviour over the active range of u a . It is noteworthy that, near the onset of the oscil- 
lations, the concentration of z (the flammable volatiles) and the temperature can rise 
dramatically at a relatively low ambient temperature. The amplitude and temperature 
range of the oscillations is, of course, also influenced by the rate of non-reactive heat 
removal i and the rate at which vapour phase material is removed from the system. 
However, it is not difficult to imagine a situation of restricted heat and mass trans- 
fer where this sudden surge in the temperature and concentration of volatiles leads to 
spontaneous ignition of the volatiles, whence the system enters a new regime of flaming 
combustion. 
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In (c) a higher rate of reactant supply induces multiple steady states, with a jump 
occurring at a comparatively low ambient temperature. Although the upper steady 
state occurs at quite low values of u and z, it is unstable, and again u and z must jump 
to a catastrophic limit cycle. (For clarity the complete branch of limit cycles is not shown 
in (c), it terminates in a homoclinic orbit.) The amplitude of the limit cycles becomes 
unrealistically high, due to the assumption of an infinite pool of substrate that is built 
into this simple model. Nevertheless, the qualitative effect is similar to that in (b): the 
temperature and volatile concentration can surge uncontrollably when the competing 
reaction is initially encouraged. 



5 Discussion and conclusion 

To a large extent, the development of fire-retarding treatments of and strategies for cel- 
lulosic materials has been based on the large body of experimental evidence, some of 
which was cited above, for the competitive nature of thermal decomposition. If char is 
formed at the expense of volatile fuel, then intuition suggests that we should expect to 
see a positive correlation between char yield and fire resistance. Many different addi- 
tives to cellulosic materials have been found to enhance char-formation, e.g., chromated 



copper arsenate treatments (CCA) |j37| , metal carboxylates [ 38fl , sodium hydroxide [32 



potassium chloride phosphates PD| , ammonium salts JIT], and those studied in |I2 



and |18[ . Some of these are used commercially. 

However, from the point of view of promoting fire safety in real situations, the idea that 
better charring properties equate to better flame resistance is too simplistic, because of 
the nonlinear feedback effects on the competition between volatile and char formation 
described in section 4. The results of this study indicate that char-forming treatments 
should be used circumspectly, with proper consideration of the particular situation of 
potential thermal decomposition — in particular, whether heat and volatiles are removed 
efficiently enough to prevent the thermal feedback that promotes volatilization. Even 
without taking into account the detailed role of the hydrolysis chemistry described in sec- 
tion ^, it is clear from thermokinetic considerations alone that there is a delicate balance 
between the amount of char and the degree of fire resistance, because the temperature 
- and therefore the rate of volatile fuel formation — can increase dramatically as the 



exothermal charring reactions occur. In [32 it was reported that pyrolysis of cellulose be- 



came strongly exothermal when the material was treated with large amounts of charring 
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retardant. Another mechanism by which extensive charring could become dangerous was 



described in j|3|: since char has much higher permeability than the unreacted substrate, 
flammable volatiles would flow more rapidly through hot char as it forms. The volatiles 
may either undergo secondary exothermal charring in the existing char, heating up the 
system further, or ignite on contact with oxygen. 

As pointed out in ||, many other factors such as rate of external heating, heat and 
mass diffusion, secondary chemistry, and the immediate environment can influence the 
decomposition process. None of these is built in to the simplified, homogeneous model 
in section f|, which is a crude and much-reduced simplification of a complex thermoki- 
netic system. Nevertheless, it does reflect the reality of the effects of cellulose thermal 
decomposition (there is nothing fine or subtle about a catastrophic spontaneous fire in a 
haystack, either), and it is consistent with the known important features of the process 
- namely, the competitive pathways, the involvement of water, and the thermokinetics. 
It is also not difficult to envisage other situations where dangerous self-heating could 
take place. It is often a statutory requirement that cellulosic materials in bedding, fur- 
nishings, and textiles in nursing homes contain fire-retardants. However, a bedding fire 
started by a cigarette is exactly the kind of situation we have in mind — the necessary 
moisture would already be present and the heat of charring cannot easily escape, thus 
provoking the volatilization reactions. 

This study of a simple model does not pretend to be quantitatively precise. Rather, 
its intent is to elicit important qualitative information concerning the dynamical be- 
haviour of a cellulosic system under thermal stress, and highlight some thermokinetic 
effects that could be a consideration in the design of targeted fire-retarding strategies for 
cellulose. Further modelling studies will examine the thermokinetics in the presence of 
fire-retardants, and simulate cellulose decomposition behaviour in an oxidizing environ- 
ment, looking at critical conditions for ignition of the volatile products and the question 
of whether decompostion is coupled to the flame chemistry as well as thermal feedback. 
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Appendix 

Eqs have the following dimensioned form, with the symbol 
and dimensionless groups as defined in the accompanying table. 



VC av — = V{-I\H l )k l e- E ^ RT c w + V(-AH 2 )k 2 e~ E ^ RT + L(T a -T) . (g) 

Symbol Definition Units 

V volume of relevant zone m 3 

c w concentration of water in zone mol/m 3 

c z concentration of volatiles in zone mol/m 3 

t time s 

T, T a system temperature, ambient temperature K 

k\ pseudo first-order dehydration rate constant s _1 

k 2 pseudo zeroth-order volatilization rate constant mol/(m 3 s) 

Ei dehydration reaction activation energy kJ/mol 

E 2 volatilization reaction activation energy kJ/mol 

R w constant rate of supply of water from reaction |l| mol/s 

F rate of outflow of vapour-phase species m 3 /s 

C av average volumetric specific heat of reacting system J/(m 3 K) 

AH reaction enthalpy kJ/mol 

L combined heat transfer coefficient J/(s K) 

R gas constant J/(molK) 

c re f a reference concentration mol/m 3 

w c w /c re{ 

u RT/Ei 

t tki 

v k 2 /(k 1 c ie {) 

H E 2 /E 1 

a AF 2 /Ai?i 

C C av E\j (c rc fR(—AHi)) 

r R w /(Vkic rel 

f F/iyh) 

i LEt/iVCrefkiRi-AH!)) 
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The carbonium ion (A) may form a levoglucosan end (B) via intramolecu- 



lar nucleophilic attack, or a reducing end (C) when intercepted by a water 



molecule. In both cases a non-reducing end (D) is also formed. 



2 How chemical and thermal feedback may be involved in cellulose decom- 


position. A, Z, W, Y, G, 


and A stand for cellulose, volatiles such as 


levoglucosan, water, dehydrated product or char, gases such as CO and 


C0 2 , and heat, respectively 





3 Steady states of Eqs. 4H 6| as a function of u a , the ambient temperature, 
for 3 values of r, the supply rate of water, (a) r = 0.9 x 10~ 4 , (b) r = 
1.1 x HT 4 , (c) r = 3.0 x 10" 4 . Other parameters: fi = 1.3, a = 0.3, 
/ = 0.0001, t = 0.004, C = 1, v = 1. Solid lines represent stable steady 
states, dashed lines are unstable steady states, and the dotted lines are 
the amplitude maxima of the limit cycles 




Figure 1: The carbonium ion (A) may form a levoglucosan end (B) via intramolecular 
nucleophilic attack, or a reducing end (C) when intercepted by a water molecule. In 
both cases a non-reducing end (D) is also formed. 
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Figure 2: How chemical and thermal feedback may be involved in cellulose decomposi- 
tion. A, Z, W, Y, G, and A stand for cellulose, volatiles such as levoglucosan, water, 
dehydrated product or char, gases such as CO and CO2, and heat, respectively. 
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Figure 3: Steady states of Eqs. [§-[! as a function of u a , the ambient temperature, for 
3 values of r, the supply rate of water, (a) r = 0.9 x 1CT 4 , (b) r = 1.1 x 10~ 4 , (c) 
r = 3.0 x 10~ 4 . Other parameters: fj, = 1.3, a = 0.3, / = 0.0001, i = 0.004, = 1, 
v = 1. Solid lines represent stable steady states, dashed lines are unstable steady states, 
and the dotted lines are the amplitude maxima of the limit cycles. 



